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Summary 

The  objective  of  this  research  was  to  develop  vibration  suppression  strategies  for  space  structures 
that  contain  frames  and  trusses  as  well  as  soft  elements.  Based  on  our  success  in  earlier  work,  our 
initial  focus  was  on  cable  control  of  frames  and  trusses  and  theoretical  and  experimental  results 
showed  good  vibration  attenuation  through  a  combination  of  two  distinct  effects:  a  parametric 
effect  that  changes  structural  stiffness  and  a  direct  effect  that  arises  due  to  external  forcing.  Cable 
actuators  do  not  lend  themselves  well  to  control  of  soft  structures  and  therefore  the  latter  part 
of  our  research  focused  on  a  new  method  of  vibration  control,  where  energy  redistribution  and 
dissipation  would  be  achieved  through  sequential  application  and  release  of  constraints.  It  was 
shown  that  application  of  constraint  typically  removes  kinetic  energy  from  the  system  and  release 
of  constraint  resets  the  system  for  a  new  cycle  of  constraint  application  and  removal.  Conditions 
that  lead  to  a  net  loss  of  kinetic  energy  per  cycle  and  bounds  on  the  amount  of  energy  removed  were 
examined.  The  experimental  component  of  this  work  included  energy  redistribution  in  a  beam  with 
electromagnetic  brakes  that  were  periodically  turned  on  and  off;  and  theoretical  extension  of  the 
work  to  soft  structures  first  looked  at  the  problem  of  passive  vibration  damping  in  a  string  using  a 
obstacle  at  its  boundary.  For  the  string  vibration  problem,  damping  can  be  increased  significantly 
by  using  a  scabbard-like  actuator  at  one  boundary  but  this  active  method  of  constraint  application 
and  removal  can  also  add  energy  to  the  system. 

1.  Vibration  Control  of  Frames  and  Trusses  Using  Cable  Actuators: 

The  numerical  model  of  a  structure  with  cable  actuators  is  developed  using  finite  elements  (Issa, 
et  al.}  2010),  and  it  takes  the  form 


MX  +  (K-  T  Kg)X  —  T  F  (1) 

where  M  and  K  are  the  mass  and  stiffness  matrices  of  the  system,  F  is  the  external  force  vector 
corresponding  to  unit  tension  in  the  cable,  X  is  the  vector  of  nodal  displacements  and  rotations, 
and  Kg  is  the  geometric  stiffness  matrix  corresponding  to  unit  cable  tension.  The  total  geometric 
stiffness  effect  from  cable  tension  is  of  the  form  TKg  since  axial  loading  in  the  frame  elements  and 
the  stiffness  effects  due  to  the  cable  elements  are  both  proportional  to  T.  It  is  useful  to  express 
the  equations  of  motion  in  terms  of  modal  amplitudes  r).  The  model  is  then  truncated  to  n  modes 
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Figure  1.  Feedback  control  design  for  frame  structures 


which  results  in  the  following  equations 

7 j  +  r tf  +  (fi  —  T  kg)rj  =  T  f  (2) 

where  kg  is  the  n  x  n  modal  geometric  stiffness  matrix  and  fi  =  diag(u;?)  is  the  diagonal  matrix  of 
zero-tension  natural  frequencies.  To  account  for  damping,  modal  damping  Trj  is  introduced  where 
T  =  diag(2£u;t)  is  the  diagonal  modal  damping  matrix  and  £t’s  are  the  modal  damping  ratios.  In 
Eq.(2),  /  is  the  modal  forcing  vector.  The  state  space  form  of  Eq.(2)  is  given  by 

X\  =  X2 

±2  =  -ftxi  -  VX2  +  (kgX  1  +  f)u 

where  X\,X2  €  Rn  are  the  modal  amplitude  and  velocity  vectors,  and  tension  in  the  cable  is  taken 
to  be  the  input,  z.e.,  u  =  T.  We  assume  the  control  input  to  be  comprised  of  a  bias  tension  and 
an  active  term  u  that  vanishes  as  the  vibrations  decay  to  zero,  i.e.,  u  =  u  -I-  T5.  The  bias  tension 
shifts  the  equilibrium  configuration  from  {x^x^)  =  (0,0)  to  (x^x^)  =  (xi,0),  where 

XI  =  [n-Tbkg]~lfTb  (4) 

where  [H  —  T5  /cs] ~ 1  exists  since  T5  is  less  than  the  critical  buckling  load.  To  shift  the  equilibrium 
configuration  back  to  the  origin,  we  define  the  following  transformation  of  coordinates  and  input 


(3) 


9i  =  (xi-xi),  q2  =  x2,  u  =u-Tb 
The  state  space  description  in  Eq.(3)  can  now  be  written  as  follows 

Qi  =  92 

q2  =  -tlqi-rq2+[kg(q1  +  xi)  + f]u  +  Tbkgqi 


(5) 

(6) 


For  the  system  described  by  Eq.(6),  the  stability  properties  of  this  equilibrium  point  (qf,  q^)  =  (0, 0) 
is  described  with  the  help  of  the  following  theorem: 

Theorem  1: 

Consider  the  frame  structure  modeled  by  Eq.(6)  with  bias  tension  Tb,  Tb  >  umin  >  0.  The  control 
design  in  Fig.l,  with  the  output  defined  by  the  relation 

y  =  -(q2  +qf  Q )  [k9(qi  +  xi)  +  /]  (7) 

renders  the  equilibrium  (qf,q2)  =  (0,0)  asymptotically  stable  for 
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Figure  2.  A  frame  structure 


1.  any  choice  of  Q  =  QT  €  Rnxn  which  guarantees  the  positive-definiteness  of  the  two  matrices 


A  = 


ft 

Q 


?)■ 


R-l  2 

B-‘  irQ  r 


l«r\ 


(8) 


2.  and  the  bias  tension  T&  satisfying 


0  <  Umiri  <  <  0 


^min(B) 


A  max(W) 

for  A max(W)  >  0,  where  A max(W)  is  the  largest  eigenvalue  of  matrix  W,  defined  as 

W  =  [  9*^ 

\k9  0 


(9) 


and  Amin(i?)  is  the  smallest  eigenvalue  of  B  and  6  €  (0, 1)  is  any  constant. 

For  numerical  simulations,  we  consider  the  frame  structure  model  depicted  in  Fig.2,  formed  by 
welding  together  aluminum  bars  of  circular  cross-section  of  diameter  6  mm  and  6  =  1.0  mm. 

The  control  model  utilized  the  first  two  modes  of  the  frame,  which  correspond  to  first  mode  bending 
in  the  xz  plane  and  simple  torsion  about  the  x  axis.  The  natural  frequencies  of  these  modes  were 
computed  using  FEA  as  u>\  =  7.92  rad/ sec  and  u 2  =  15.20  rad/ sec,  respectively.  The  damping 
ratios  of  these  modes  were  taken  to  be  £1  =  £2  =  0.003.  To  ensure  positive-definiteness  of  matrices 
A  and  B  in  Eq.(8),  Q  was  chosen  as  Q  =  0.0157T.  In  accordance  with  Eq.(9),  the  bias  tension  and 
y>min  were  chosen  as  =  2  umin  =  0.18  N.  For  the  cable  placement  shown  in  Fig.2,  the  geometric 
stiffness  matrix,  the  modal  forcing  vector  and  the  equilibrium  configuration  were  computed  (in  SI 
units)  and  found  to  be  as  follows 


f  0.326  0.319  \ 

V  0.319  —0.328  )  ’  1 


x  10~5 


Simulation  results  are  shown  in  Fig.3  for  initial  conditions  <?i(0)  =  (0,0)  and  92(0)  =  (0.12,0.12) 
in  SI  units,  and  control  system  parameters  of  r  =  1/15  sec,  e  =  1  x  10“7,  and  umax  =  30  N.  It  is 
clear  from  the  plots  that  the  vibration  of  the  modeled  modes  are  very  effectively  controlled. 

Experiments  were  carried  out  in  the  laboratory  with  a  frame  similar  to  the  one  in  Fig.2.  The 
experimental  results  can  be  found  in  the  paper  by  Issa,  et  al.  (2010). 
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Figure  3.  Plot  of  modal  amplitudes  q\  =  (gi,G2)T  and  control  effort 

2.  Energy  Dissipation  Through  Application  and  Removal  of  Constraints: 

Consider  the  iV-dof  linear  dynamical  system 

MX  +  KX  =  0  (10) 

where  M  and  K  are  the  iV-dimensional  mass  and  stiffness  matrices  and  X  =  [xi,X2,  •  •  •  ,xat]T  is 
the  vector  of  independent  generalized  coordinates.  Upon  application  of  a  holonomic  constraint,  the 
dynamics  of  the  system  takes  the  form 


MY  +  KY  =  0  (11) 

where  M  and  K  are  the  ( N-. 1  )-dimensional  mass  and  stiffness  matrices  and  Y  =  [yi ,  J/2 >  *  •  •  ,  y(AT-i)] T 
is  the  vector  of  independent  generalized  coordinates  of  the  constrained  system. 

We  denote  the  unconstrained  system  as  state  a  and  the  constrained  system  as  state  (3.  The 
transition  from  state  a  to  state  (3  occurs  over  the  brief  interval  of  time  when  the  constraint  is 
applied.  If  \t~p,  t+p]  denotes  this  interval,  the  state  transition  can  be  described  by  the  relations 


= 

(12) 

a 

=  MX(t~p)  +  Iq-,0 

(13) 

=  To*  *(«£,) 

(14) 

*■(*&) 

=  Ta0X(t+0) 

(15) 

where  Ia->/3  is  the  iV-dimensional  impulse  of  the  generalized  forces  and  Tap  is  a  transformation 
matrix.  The  transition  from  state  (3  to  state  a  occurs  over  a  brief  interval  of  time  when  the 
constraint  is  removed.  If  [tja,  denotes  this  transition  interval,  the  effect  of  the  transition  can 
be  described  by  the  relations 


=  nw 

(16) 

ntt) 

=  nt-f,,) 

(17) 

*(#«) 

=  TgaY(t%,) 

(18) 

x{t%) 

=  TtaY(t+J 

(19) 

where  Tpa  is  a  constant  and  unique  transformation  matrix.  The  change  in  kinetic  energy  over  one 
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cycle  of  constraint  application  and  removal  is 

AE  =  \xT{tt0)MX{tt0)  ~  \xT{t-mfi)MX{t~&)  (20) 

and  through  simplifications,  it  can  be  shown  to  be  equal  to 

AE  =  \xT{t~ <  0  (21) 


In  modal  coordinates,  the  change  in  kinetic  energy  can  be  expressed  as  follows 


AE  =  \XT{tZp)M[x{tt0)-X{t-a(i)] 

=  |/iT(f^)^TM^I>(f:/3)-$A(W] 

=  ~  5A7'(^q/3)-^A(^q^)  —  0 

(22) 

where 

a  =  (iN-  rTr) 

(23) 

r  = 

(24) 

$  =  [<t>u<h,- ■  ■  ,M  e  Rn*n  and  #=[&,&,•••  ,V_d]  €  RN*(N~') 

are  modal  matrices  in 

states  a  and  /?,  respectively,  and  Ijv  is  the  iV-dimensional  identity  matrix. 


To  understand  the  limitations  of  energy  dissipation  through  application  of  constraints,  consider 
the  matrix  TTr,  where  T  is  defined  by  Eq.(24).  It  can  be  shown  that  TtT  is  idempotent  and  hence 
A  =  (J/v  —  rTr)  is  idempotent.  The  trace  of  A  can  be  computed  as 

trace[A]  =  trace[J/v]  —  tracepT^T]  =  N  —  tracepTT71]  =  N  —  trace  =  1  (25) 

From  Eq.  (25)  we  deduce  that  A  has  one  unit  eigenvalue  and  other  eigenvalues  are  all  zero.  If  v  is 
the  normalized  eigenvector  of  A  corresponding  to  its  unit  eigenvalue,  then  eigen  decomposition  of 
A  gives 

A  =  vvT  (26) 

Substitution  of  Eq.(26)  into  Eq.(22)  gives  the  expression  for  the  loss  in  kinetic  energy 

=  ~\iiT {t~0)wT [«r/i(Q)]  <  0  (27) 

It  is  clear  from  Eq.(27)  that  the  energy  of  the  system  will  not  be  dissipated  if  any  one  of  the 
following  conditions  hold: 

1.  =  0:  the  constraint  is  applied  when  the  system  has  zero  kinetic  energy. 

2.  vT  =  0,  Vi  7^  0,  i  =  1,2,  •••  ,AT:  the  constraint  is  applied  when  the  modal  velocity 

vector  is  normal  to  the  eigenvector  of  A  corresponding  to  the  unity  eigenvalue. 

3.  Vi  =  0,  i  G  Sr  ==  {/ci,  •  •  •  ,  /cr},  Aj(^)  —  V?  ^  £r:  the  kinetic  energy  of  the  system  lies 

in  specific  modes  that  correspond  to  zero  entries  of  v. 

The  first  and  second  conditions  can  be  avoided  through  a  proper  choice  of  time  when  the  constraint 
is  applied.  The  third  condition  corresponds  to  a  necessary  condition  for  energy  entrapment  and  can 
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occur,  for  example,  when  the  constraint  is  applied  at  a  node  of  a  system.  The  sufficient  condition 
for  energy  entrapment  is  stated  next  with  the  help  of  the  following  theorem. 

Theorem  1  The  energy  of  a  linear  system  will  remain  trapped  in  r  specific  modes  if  T  contains 
an  orthonormal  sub-matrix  of  dimension  r. 

For  numerical  simulations,  consider  the  system  in  Fig.4  (Issa,  et  al.,  2009).  For  this  system,  T, 
A  and  v  were  computed  as  follows: 

/  0.9599  0.0  \  /  0.0786  0.0  0.2691  \  /  0.2804 

rr=  0.0  1.0  ,  A=  0.0  0.0  0.0  ,  v=i  0.0 

\  -0.2804  0.0  )  \  0.2691  0.0  0.9214  )  \  0.9599 


-^r  rotational  spring  -(o)-  free  pin  joint  locked  pin  joint 


Figure  4.  A  linear  system  with  3DOF  in  state  a  and  2DOF  in  state 


From  the  entries  of  T  in  Eq.(28)  it  is  clear  that  the  condition  in  Theorem  1  is  satisfied  with  r  =  1. 
Since  =  1-0,  the  energy  of  the  system  is  trapped  in  the  second  mode  of  state  a,  which  is 
also  the  second  mode  of  state  /?.  We  can  also  verify  that  one  element  of  u,  namely  t>2,  is  equal  to 
zero.  Energy  entrapment  can  be  verified  from  the  simulation  results  presented  in  Fig.5,  the  initial 
conditions  for  which  were  chosen  as 

(xi,X2,X3,ii,i2)i3)  =  (0.006,0.016,0.002,0.00,0.00,0.00) 

in  SI  units.  In  this  simulation,  the  system  was  switched  from  state  a  to  state  f3  at  the  earliest 
opportunity  after  0.2  seconds  when  the  two  middle  bars  are  aligned,  i.e.,  when 

g(X)  =  xi  -  2x2  4-  x3  =  0 

and  switched  back  to  state  a  after  0.2  seconds  in  state  /?.  The  modal  amplitudes  in  state  a,  namely, 
Mi,  M2?  M3?  are  shown  in  Fig.5  for  the  intervals  of  time  when  the  system  is  in  state  a.  It  is  clear  from 
these  plots  that  the  amplitudes  of  the  first  and  third  modes  decay  to  zero  whereas  the  amplitude 
of  the  second  mode  remains  constant.  The  plot  of  the  energy  confirms  that  some  energy  of  the 
system  gets  trapped  in  the  second  mode.  This  is  not  surprising  since  the  location  of  the  brake  at 
the  node  of  the  second  mode  renders  it  ineffective. 

3.  Modal  Disparity  and  its  Experimental  Verification: 

For  experimental  validation  of  the  work  presented  in  the  previous  section,  we  investigated  the 
free  vibration  of  a  beam  with  a  mid-span  hinge,  as  shown  in  Fig.6.  Assuming  Euler-Bernoulh  beam 
theory,  the  equation  of  motion  of  the  beam  in  the  x-y  plane  can  be  written  as  follows 

El y""  +  my  =  0  if  x  €  (0, L/2)  or  x  G  (L/ 2,  L)  (29) 
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Figure  5.  Plot  of  modal  amplitudes  in  state  a  and  total  energy  of  the  system 


and 


2/(0,  t)  =y(L,t )  =  0 
y(0,t)  =  y'(L,t)  =  0 
y((L/2)-,t)  =  y((L/2)+,t) 
y'"((L/2)-,t)  =  y"\(L/2)+,t) 


y"((L/2)-,t)  =  y"((L/2)+,t)  =  0 


in  stiffness  state  a:  hinge  released 


(30) 


y'((L/2)-,t)  =  y((L/2)+,t) 
y"  {{L/2)~  ,t)  =  y"((L/2)+,t] ) 


in  stiffness  state  /? :  hinge  locked 


(31) 


F 
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Figure  6.  A  flexible  clamped-clamped  beam  hinged  in  the  middle 


In  Eq.(29),  E  is  the  modulus  of  elasticity,  I  is  area  moment  of  inertia,  and  m  is  the  mass  per 
unit  length  of  the  beam.  The  beam  is  modeled  using  N  standard  (cubic)  finite  elements,  with  two 
degrees  of  freedom  per  node  (translation  and  rotation  about  the  2  axis).  To  facilitate  the  modeling 
of  the  hinge  in  its  two  states,  the  node  at  x  =  L/2  has  two  rotational  degrees  of  freedom,  6l  and  0r , 
corresponding  to  ?/(L/2~,  t)  and  y  (L/ 2+,  £),  respectively.  In  stiffness  state  a,  the  hinge  is  assumed 
to  be  free  and  6l  and  6r  are  independent.  However,  in  stiffness  state  /?,  the  hinge  is  locked  and 
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the  constraint  9l  =  9r  has  to  be  enforced.  This  is  accomplished  by  adding  a  penalty  of  magnitude 
^ kr  ( 9l  —  6r  )2  to  the  strain  energy,  that  is,  by  adding  a  rotational  stiffness 


Kr  =  kr 


1 

-1 


(32) 


to  the  global  stiffness  matrix.  The  parameter  kr  is  set  to  zero  in  the  stiffness  state  a  and  to  a  large, 
positive  value  in  stiffness  state  /?.  With  the  additional  degree  of  freedom  (in  rotation)  for  the  node 
at  x  =  L/2,  the  finite  element  model  has  ( 2N  —  1)  degrees  of  freedom. 

A  lumped  mass  M  is  added  at  the  central  node  to  account  for  the  mass  of  the  hinge  and  the 
electromagnetic  brake.  The  transition  from  stiffness  state  a  to  stiffness  state  (3  is  achieved  by 
activating  an  electromagnetic  brake.  The  activation  of  the  brake  occurs  over  a  brief  interval  of  time 
and  results  in  the  application  of  an  action-reaction  pair  of  impulsive  moments  to  the  middle  node. 
If  t  G  [tapytap]  denotes  the  activation  time,  the  effect  of  the  impulsive  moments  can  be  described 
by  the  relations 


Y(tZfi)  =  Y(fo) 

(33 

MY(t~0)  +  IQ^  =  MY(rt +,) 

where  M  is  the  mass  matrix  and  Y  is  the  vector  of  nodal  degrees  of  freedom.  Ia-*p  is  the  impulse 
vector  with  nonzero  entries  corresponding  to  the  coordinates  9l  and  9r  and  has  the  form 


Ja-tf  =  [<>,•••  ,0]T, 


(34) 


where  C  is  the  impulse  and  r  is  the  impulsive  moment. 

The  transition  from  stiffness  state  (3  to  stiffness  state  a  is  achieved  by  releasing  the  brake.  If 
t  G  [tpQl  tpa]  denotes  the  brief  time  interval  over  which  the  brake  is  released,  the  degrees  of  freedom 
and  their  velocities  just  prior  to  and  after  release  of  the  brake  are  the  same,  i.e. 


Y(tpQ)  =  Y{t+0a),  Y{t-0Q)  =  Y{t+a)  (35) 

With  this  notation,  the  stiffness  parameter  kr  can  be  defined  as  follows 

^  _  f  0  stiffness  state  a 

r  \  koo  stiffness  state  /? 

where  k^  is  some  large  positive  number  chosen  to  enforce  the  constraint  9l  =  9r. 

Consider  the  clamped-clamped  beam  in  Fig.6  with  the  following  properties: 

Table  1.  Material  and  geometric  properties  of  clamped-clamped  beam  in  Fig.6 


Material 

Aluminum 

Young’s  Modulus 

71  GPa 

Density 

2710  Kj/rr? 

Dimensions 

2.0  x  0.05  x  0.0023  m 

Hinge  mass 

0.182  Kg 

The  natural  frequencies  of  the  first  four  modes  of  the  beam  in  the  two  stiffness  states  axe  shown  in 
Table  2.  Since  the  hinge  is  located  at  mid-span,  the  natural  frequencies  of  even  numbered  modes  in 
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the  two  stiffness  states  are  identical.  This  is  true  because  even-numbered  modes  have  zero  curvature 
at  mid-span  and  are  not  affected  by  the  state  of  the  hinge,  i.e.  locked  or  released.  The  first  four 
mode  shapes  of  the  beam  in  the  two  stiffness  states  are  provided  in  Fig.7  for  reference. 

Table  2.  Natural  frequencies  of  the  finite  element  model  of  the  beam  in  the  two  stiffness  states 


uai,  vpi  (rad/sec) 

mode  number,  i 

i  =  1 

i  =  2 

i  =  3 

i  —  4 

stiffness  state 

a 

1.29 

8.34 

9.52 

27.0 

0 

2.30 

8.34 

14.68 

27.0 

Using  modal  truncation,  the  matrix  measure  of  modal  disparity  was  computed  using  the  first  four 


modes  as  follows 


0.980 

0.000 

0.153 

0.000  ^ 

0.000 

1.000 

0.000 

0.000 

0.139 

0.000 

0.949 

0.000 

0.000 

0.000 

0.000 

1.000  j 

(36) 


The  second  and  fourth  rows  and  columns  of  the  matrix  in  Eq.(36)  maintain  the  identity  structure. 
This  is  indicative  of  the  fact  that  even-numbered  modes  in  the  two  stiffness  states  are  identical. 
The  nonzero  elements  in  the  off-diagnonal  entries  of  the  matrix  are  indicative  of  the  presence  of 
modal  disparity  between  the  two  stiffness  states  and  indicate  how  modal  energy  will  be  redistributed 
between  odd-numbered  modes  in  these  two  states. 

Consider  the  scenario  where  the  beam  is  initially  in  stiffness  state  a  and  vibrating  purely  in  the 
third  mode  with  an  amplitude  A.  The  total  energy  of  the  beam  is  equal  to  Ea  =  0.5A2  It  is 
assumed  that  there  is  no  damping  in  the  system  and  that  the  brake  is  activated  when  the  beam 
passes  through  its  neutral  position.  The  modal  velocities  in  stiffness  states  ft  can  be  computed  as 


( 

^(C/s) 

*><(*Jj) 

\ 

.980 

0.000 

0.153 

.000 

1.000 

0.000 

.139 

0.000 

0.949 

.000 

0.000 

0.000 

0.000  \ 

'  0  ' 

0.000 

0 

0.000 

Auja2 

1.000  ) 

0 

(37) 
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Clearly,  the  energy  of  the  beam  is  redistributed  in  modes  1  and  3  in  stiffness  state  (3.  The  amount 
of  energy  in  these  modes  are 

Efii  =  |  0.1532j42  cj23  =  0.1532  Ea,  Efr  =  \  0-9492.42  cj23  =  0.9492  Ea  (38) 

It  can  be  easily  shown  that  the  amplitudes  of  these  modes  are  OAbZA^^/upi)  =  0.633A  and 
0.§4§A(uaz/u)pz)  =  0.615A,  respectively.  These  results  have  been  validated  through  experiments. 
A  description  of  the  experimental  setup  can  be  found  in  the  paper  by  Issa,  et  al.,  (2008). 

4,  Vibration  of  a  String  Wrapping/Unwrapping  an  Obstacle  in  its  Boundary: 

Consider  a  string  vibrating  against  an  obstacle  placed  at  one  of  its  boundaries,  as  shown  in 
Fig.8.  We  investigate  energy  dissipation  in  the  string  under  the  following  assumptions: 


y 


Figure  8.  A  string  vibrating  against  an  obstacle. 


Al.  The  obstacle  is  rigid  and  has  the  following  geometry 

y  =  f(x),  y(  0)  =  0, 


A2.  The  string  is  homogenous  and  has  a  constant  mass  per  unit  length  denoted  by  p.  The  tension 
in  the  string  is  equal  to  T  and  remains  constant  at  all  times.  The  string  undergoes  transverse 
vibration  in  the  xy  plane  and  is  not  affected  by  gravity. 

A3.  The  amplitude  of  oscillation  of  the  string  is  small  and  therefore  the  equation  of  motion  of  the 
string  can  be  expressed  by  the  standard  relation 


where  y(x,t)  is  the  displacement  of  the  string  at  a  distance  x  from  the  origin  at  time  t. 

A4.  The  string  wraps  around  the  obstacle  during  vibration.  Over  each  time  step  during  wrapping, 
a  small  element  of  the  string  comes  to  rest  on  the  obstacle  through  perfectly  inelastic  collisions. 
The  wrapping  process  continues  till  the  freely  vibrating  portion  of  the  string  has  no  more 
kinetic  energy. 

A5.  The  surface  of  the  obstacle  is  not  sticky  and  the  string  unwraps  from  the  obstacle  without 
any  loss  of  energy. 

A6.  At  the  initial  time  t  =  0,  the  string  has  no  contact  with  the  obstacle.  It  is  in  its  mean  position 
with  zero  potential  energy  and  kinetic  energy  equal  to 


dy 

dx 


=  0 


T.=0 


(39) 
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A7.  The  string  continues  to  vibrate  in  the  mode  in  which  it  started  its  vibration  at  the  initial 
time.  This  implies  that  each  point  of  the  string,  not  in  contact  with  the  obstacle,  has  the 
same  frequency  of  vibration  at  any  instant  of  time,  and  the  number  of  nodes  in  the  vibrating 
string  remains  constant. 

A8.  The  string  has  no  internal  damping,  i.e .,  the  energy  of  the  string  will  remain  conserved  during 
free  vibration. 


A  general  solution  to  the  partial  differential  equation  in  Eq.(39)  can  be  written  as 

y(x,  t)  =  A  (sin  Ax  —  tan  XI  cos  Ax)  sin  cut 


(41) 


We  now  consider  the  boundary  conditions  at  the  contact  break  point.  From  Fig.10  we  have 


f(x)  =  y(x,  t)  =>  f(x)  =  A  (sin  Ax  —  tan  XI  cos  Ax)  sin  ut  (42) 

Also,  the  string  is  tangential  to  the  obstacle  at  the  contact  break  point  x  =  x,  i.e., 
dy 

f'(x)  —  —  (x,  t)  =»  /'(x)  =  AA  (cos  Ax  +  tan  A l  sin  Ax)  sin  ujt  (43) 

ox 

From  Eqs.(42)  and  (43)  we  get 

tanA(i-x)  =  (44) 

which  indicates  that  A  can  be  computed  from  the  value  of  x.  The  solution  of  Eq.(44)  is  however 
not  unique  -  each  non-trivial  value  of  A  corresponds  to  a  mode  of  vibration  of  the  string. 

Let  the  total  energy  of  the  string  at  any  time  t  be  denoted  by  E.  Then, 


E  =  Epe  +  Eke  =  Eg*  +  ET*  +  Eke  (45) 

where  E°^s  is  the  potential  energy  of  the  string  wrapped  around  the  obstacle,  is  the  potential 
energy  of  the  freely  vibrating  string,  jEpe  is  the  total  potential  energy,  and  E^e  is  the  kinetic  energy 
of  the  string.  The  total  potential  energy  of  the  string  is  computed  as  the  product  of  the  tension  T 
and  elongation  of  the  string 


Epe=T  J  dl  =  T  J  (ds  —  dx)  =  T  J  (y/dx2  +  dy2  -  dx)  =  T  t (</!  +  {dy/dx)*  -  l)dx 


and  £peS  and  can  be  written  as 


£pes  =  T  /Vl  +  {dy/dxY  -  1  )dx,  Eg>  =  T  f\y/l  +  (dy/dx)2  -  1  )dx 
Jo  Jx 

Since  the  string  conforms  to  the  shape  of  the  obstacle,  (dy/dx)  =  f'(x)  for  x  £  [0,  x].  By  expressing 
(dy/dx)  =  y'(x,  t)  for  x  G  [x,Z]  and  simplifying,  we  get 


£p°eS  =  [^1 +  [/'(*)] 

f  =  Tf[/: 


1 


dx 


(46) 


77<vib 

■^pe 


l  +  [j/(x,t)]  -1 


dx  =  ^TAsec2  Al|2A(l  —  x)  +  sin[2A(Z  —  5)]|  [5(x)]2(47) 
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Figure  9.  Exponential  decay  in  the  energy  of  a  string  wrapping  and  unwrapping  around  an  obstacle.  The 
plots  show  energy  decay  for  single-mode  vibration  in  the  first  four  modes  with  Eq  =  0.5  J. 

where  g(x)  is  defined  in  Alsahlani  and  Mukherjee  (2010(a)).  The  kinetic  energy  of  the  vibrating 
string  can  be  written  and  simplified  as  follows 

£ke  =  \  J£  p[if{x,t)}2dx  =  -£)-sin[2A(Z  [p(^)]2}  (48) 

From  the  above  equations  it  is  easy  to  verify  that  the  energy  expression  has  the  form 

E  =  h(x,A)  (49) 

For  a  configuration  in  which  the  string  is  wrapped  around  the  obstacle,  the  complete  solution  can 
be  determined  from  the  values  of  x  and  E  using  the  four-step  algorithm  described  in  Alsahlani  and 
Mukherjee  (2010(a)). 

For  numerical  simulations,  consider  a  string  with 

T  =  IN,  p  =  0.025  kg/m,  i  =  4m  (50) 

The  obstacle  is  assumed  to  be  a  circle  of  radius  R  and  center  coordinates  (x,y)  =  (0 ,i2),  t.e., 

y  =  f(x)  =  R~  \J R2  —  x2,  0  <  x  <  R  (51) 

Figure  9  plots  the  decay  in  energy  as  a  function  of  time  for  vibration  in  the  first  four  modes  with 
Eq  =  0.5  J.  The  following  observations  can  be  made: 

For  any  mode  of  oscillation,  it  can  be  seen  that  the  percentage  energy  loss  is  higher  for  higher 
values  of  Eq.  This  is  not  surprising  since  higher  values  of  Eq  results  in  higher  kinetic  energy  and 
greater  length  of  wrapping,  as  evident  from  the  values  of  x k  in  Table  3,  and  consequently  more 
energy  loss  through  inelastic  collision. 

The  percentage  energy  loss  is  higher  for  higher  modes  of  oscillation  for  the  same  value  of  Eq. 
This  is  true  for  the  same  number  of  cycles  as  well  as  for  the  same  length  of  time  and  is  due  to  the 
fact  that  the  velocities  of  the  string  associated  with  higher  frequencies  are  higher  in  higher  modes, 
and  as  a  consequence  the  loss  upon  impact  is  higher.  The  value  of  Xk  is  less  for  the  higher  modes 
but  this  does  not  have  a  significant  effect  on  percentage  of  energy  loss. 
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5.  Vibration  Control  of  a  String  Using  a  Scabbard-Like  Actuator: 

Consider  the  vibrating  string  in  Fig.  10(a),  that  passes  through  a  scabbard  located  at  its  left 
boundary.  At  time  t  =  tc,  the  scabbard  is  moved  instantaneously  to  the  right  by  distance  xo  along 
the  mean  position  of  the  string.  This  is  shown  in  Fig.  10(b).  At  some  future  time  t  =  trytT>  tc,  the 
scabbard  is  moved  back  to  its  original  position.  To  investigate  the  effect  of  application  and  removal 
of  the  scabbard  on  the  vibration  of  the  string,  we  assume  that  the  movement  of  the  scabbard 
imposes  a  zero-displacement  constraint  over  the  length  interval  x  G  [0,xo).  The  displacement  and 
velocity  of  the  string  over  the  remaining  interval  x  G  [xo,Z]  remains  unchanged  immediately  after 
movement  of  the  scabbard.  Additionally,  at  time  t  =  tr,  tT  >  tc,  the  scabbard  is  instantaneously 
moved  back  to  its  original  position,  f.e.,  to  the  left  by  a  distance  xo  along  the  mean  position  of  the 
string.  The  displacement  and  velocity  of  the  string  over  the  interval  x  G  [xo ,Z]  remains  unchanged 
immediately  after  movement  of  the  scabbard. 


approximate  t » tc 


y  i 

A 

t\  — 

*  dcforma 
\ 

tion^- - 

/  ^expected  real  deformation  \ 

scabbard 

1 

1 

Figure  10.  A  zero- displacement  constraint  is  applied  to  the  string  at  time  t  =  tc  over  the  length  segment 
x  G  [0,  xo)  using  a  scabbard. 

We  skip  the  mathematical  development  here  which  can  be  found  in  Alsahlani  and  Mukherjee 
(2010(b))  and  present  simulation  results  for  one  cycle  of  constraint  application  and  removal. 

We  consider  a  string  of  length  Z  =  4  m,  mass  per  unit  length  p  =  0.25  kg/m,  and  tension  T  =  1 
N,  vibrating  in  its  first  mode  with  unit  amplitude  Ao  =  1  m.  At  t  =  0  the  string  is  assumed  to 
pass  through  its  mean  position,  and  the  equation  of  motion  of  the  string  is 

7T  7T 

2/o OM)  =  sin(-x)  sin(-i)  (52) 


Figure  11.  (a)  Position  of  the  string  at  different  time  instants  (b)  Percentage  change  in  energy  of  the  string 
due  to  application  of  the  constraint  for  different  values  of  Xq  and  tc. 

The  time  period  of  oscillation  of  the  string  is  4  sec  and  Fig.  11  (a)  shows  the  shape  of  the  string 
at  different  instants  of  time  over  the  interval  t  G  [1.0, 3.0]  sec.  The  string  has  maximum  potential 
energy  at  t  =  1.0  sec  and  t  =  3.0  sec,  and  maximum  kinetic  energy  at  t  =  2.0  sec. 

We  present  simulation  results  for  percentage  change  in  energy  due  to  constraint  application 
for  xq  £  {0.0Z,  0.01Z,  0.02Z,  •  ••  ,0.09Z,0.10Z}  and  tc  G  [1.7, 2.3]  sec.  The  results,  shown  in  Fig.ll(b), 
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indicate  that  the  energy  of  the  string  increases  if  the  constraint  is  applied  when  the  string  is  far 
away  from  the  mean  position,  and  decreases  if  the  constraint  is  applied  when  the  string  is  near  its 
mean  position,  irrespective  of  the  value  of  Xq.  Upon  application  of  the  constraint,  the  potential 
energy  of  the  string  increases.  This  increase  is  large  when  the  string  is  far  away  from  its  mean 
position,  and  although  the  constraint  removes  the  kinetic  energy  of  a  portion  of  the  string,  there 
is  a  net  gain  in  energy.  When  the  string  is  near  its  mean  position,  the  increase  in  energy  due 
to  change  in  potential  energy  is  small  compared  to  the  loss  of  kinetic  energy  and  as  a  result  the 
net  change  in  energy  is  negative.  When  the  string  is  at  its  mean  position,  there  is  no  change  in 
potential  energy  upon  application  of  the  constraint  and  consequently  the  energy  loss  is  maximum 
in  this  configuration.  For  the  same  value  of  tc,  a  larger  value  of  xq  results  in  a  larger  increase  in 
potential  energy  and  a  larger  decrease  in  kinetic  energy. 

We  have  investigated  the  effect  of  application  of  the  constraint  on  the  dynamics  of  the  string 
with  arbitrary  initial  conditions  such  that  sequential  application  and  removal  of  constraints  can  be 
explored  as  a  strategy  for  vibration  suppression  -  see  Alsahlani  and  Mukherjee  (2010). 
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